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Abstract - We propose the minimally nonlinear irreversible heat engine as a new general theo- 
retical model to study the efficiency at the maximum power 77* of heat engines operating between 
the hot heat reservoir at the temperature Th and the cold one at Tc (Tc < T^)- Our model is based 
on the extended Onsager relations with a new nonlinear term meaning the power dissipation. In 
this model, we show that rf is bounded from the upper side by a function of the Carnot efficiency 
rjc = 1 — Tc/Th as rf < r)c / ["2 — Vc)- We demonstrate the validity of our theory by showing that 
the low-dissipation Carnot engine can easily be described by our theory. 



Introduction. — Facing with the recent worldwide 
problems such as the global warming and the depletion of 
energy resources, we have been urged to a low-carbon sus- 
tainable society. Demands for more efficient and greener 
heat engines have rapidly been rising since heat engines 
convert heat energy into useful work by utilizing only tem- 
perature difference, which is abundant in the earth's en- 
vironment: gcothermal power generation, solar thermal 
power generation, etc. might be promising candidates, for 
example. To evaluate and control performance of heat en- 
gines, it must be important to know the upper bound of 
the energy conversion efficiency of them. The efficiency 77 
for a heat engine is defined as 77 = W/Qh, where Qh and 
W denote the heat transferred from the hot heat reservoir 
at the temperature Th and the work output, respectively. 
Defining Qc as the heat transferred from the cold heat 
reservoir at the temperature (< T/i), we can express 
W as W = Qh + Qc- Thermodynamics tells us that rj is 
bounded from the upper side as 



Tc 

V - TTT =110, 



(1) 



where rjc denotes the Carnot efficiency and the equal- 
ity holds only when the heat engine is infinitely slowly 
(quasistatically) operated to satisfy reversibility. Because 
the heat engine realizing the Carnot efficiency takes in- 
finte time to output a finite amount of work, its power 
(work output per unit time) is absolutely and thus it 
is of no practical use. Motivated by this fact, Curzon 
and Ahlborn [l] proposed a phenomenological finite-time 
Carnot cycle model and derived that the efficiency at the 



maximum power 77* of their model is given by an appealing 
expression as 



1 



VCA, 



(2) 



which reminds us of the Carnot efficiency. Historically 
and strictly speaking, the formula 1 — y/Tc/Th itself was 
derived by others [21 [3] more previously than [1]. But 
it is usually called the Curzon-Ahlborn (CA) efficiency. 
We also call it the CA efficiency here in accord with 
the custom. The paper by Curzon and Ahlborn trig- 
gered subsequent studies on the efficiency at the maxi- 
mum power of various heat engine models [4U37j . Among 
recent studies on the CA efficiency, it is an important 
progress that Van den Brocck [12] proved that the CA 
efficiency AT/(2T) = tjca + ©(AT^) is the upper bound 
of the efficiency at the maximum power for the heat en- 
gines working in the linear response regime. Here, we 
have defined the temperature difference AT = Th — 
which is assumed to be small, and the averaged temper- 
ature T = {Th + Tc)/2, respectively. Those heat engines 
working in the linear response regime, which we call the 
linear irreversible heat engines, are described by the fol- 
lowing Onsager relations [551155] : 



Ji — LiiXi + L12X2, 
J2 = L21X1 + L22X2, 



(3) 
(4) 



where Xi = F jT^ — F/T with an external force F, Ji = i 
with the conjugate variable x oi F, X2 = l/Tc — l/Th — 
AT/T^, J2 = Qh and L^-'s are the Onsager coefficients 
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with the reciprocity L12 = L21. The dot denotes the quan- 
tity per unit time. Regarding Xi as the control parameter 
for the maximization of the power P = —Fx = —JiXiT, 
we can sec that the maximum power P* and the efficiency 
at the maximum power 77* of the hnear irreversible heat 
engines described by eqs. ([3]) and are given by 



P* = 



V = 



q^L22AT' 
AT g2 



2T 2 



(5) 
(6) 



respectively [T^], where q is called the coupling strength 
parameter and is defined as 



L12 



V L11L2 



(7) 



Since the positivity of the entropy production rate a = 
JiXi + J2X2, which is a quadratic form of the thermody- 
namic forces, restricts the Onsager coefficients L^j's to 



Lii > 0, L22 ^ 0, L11L22 ~ L12L21 > 0, 



(8) 



q should take ~l < q < 1. Thus rf in eq. ([6]) takes the 
upper bound AT/{2T) which is equal to the CA efficiency 
up to the linear order of AT, when the tight-coupling con- 
dition 



L12 



22 



1 



(9) 



holds. This condition is equivalent to saying that the two 
thermodynamic fluxes become proportional as J2 oc Ji. 
Due to the generality of the theory, the study by Van 
den Broeck renewed the interests in the CA efficiency 
and inspired the studies on the efficiency at the maxi- 
mum power of various heat engine models: it has been 
indeed shown that 77* of many heat engine models, rang- 
ing from the steady-state brownian motors [l3l[T71[T5 to 
the finite-time Carnot cycles [SOIISII, is given by eq. ^ 
in the linear response regime. Even when the system be- 
gins to enter the nonlinear response regime, it is often 
observed that rj* agrees with the CA efficiency up to the 
quadratic order of AT as rj* = AT/(2T) + AT'^/{8T^) = 
r]cA + 0{AT^) [2H171[31[3S|. This fact was firstly ob- 
served in [24] and proposed as a conjecture in [25]. Later, 
it was proved to be a precise result for the system which 
satisfies the left-right symmetry condition in addition to 
the tight-coupling condition in [27]. However, it is also 
shown, for example in [Tgi[^[MH2^l34[|35| . that ry* can 
exceed the CA efficiency in the nonlinear response regime, 
when those conditions do not hold. Therefore the CA ef- 
ficiency is no longer the upper bound of rj* for nonlinear 
irreversible heat engines and we need to construct a gen- 
eral theory to determine the upper bound of rj* for them. 

In this paper, we propose the minimally nonlinear ir- 
reversible heat engine described by the extended Onsager 
relations, where a new nonlinear term — 7^ Ji^ meaning the 



power dissipation is added to eq. (jl]) (see eq. ([T^ V The 
addition of this new nonlinear term can be seen as a natu- 
ral extension of the linear irreversible heat engine and we 
formulate 77* of such nonlinear irreversible heat engines. 
Our new formula cq. (|20p contains the coupling strength 
parameter q and a parameter jc/lh^ where 7c and 7^, de- 
note the degree of dissipation to the cold and hot heat 
reservoirs, respectively. We show that our rj* is bounded 
from the upper side by a function of the Carnot efficiency 
rjc as 1]* < ri+, where ?/+ = ijc / {2 — rjc) ■ Remarkably 
this 77+ was also mentioned in the previous studies on var- 
ious finite-time heat engine models [MUSMSSSnilll] ■ The 
generality of our theory allows us to unify these previous 
results and explain the universality of 77+. For a demon- 
stration of the validity of our theory, we show that a finite- 
time Carnot cycle model, called the low-dissipation Carnot 
engine j35j , can be described by the extended Onsager re- 
lations. 

Extended Onsager relations. Let us consider 
that a general heat engine is working between the hot and 
the cold heat reservoirs with the temperature difference 
AT = Tfi — Tc- Heat engines are generally classified into 
two types: steady-state heat engines and cyclic heat en- 
gines. Steady-state heat engines literally work in a steady 
state since an external force applied on the heat engines 
is time-independent and the hot and cold heat reservoirs 
contact with the heat engines simultaneously. Cyclic heat 
engines, in contrast, work cyclically in a time-dependent 
way such that the hot and the cold heat reservoirs con- 
tact with the heat engines alternately, not simultaneously. 
Our theory below can treat both types of heat engines in 
a unified manner. We can generally write the total en- 
tropy production rate a of the heat engine as the entropy 
increase rate of the heat reservoirs as 



Th Tc Tc [Tc Th 



(10) 



where we do not need to consider the entropy increase 
inside the heat engine since the heat engine itself is always 
in a steady state or comes back to the original state after 
one cycle. Note that the dot denotes the quantity per 
unit time for steady-state heat engines and the quantity 
divided by the one-cycle period Tcyc for cyclic heat engines. 
The power P = W is expressed as P = —Fx for steady- 
state heat engines where the time-independent external 
force F is acting on its conjugate variable x, and as P = 
W/Tcyc for cyclic heat engines. From the decomposition 
cr = JiXi + J2X2, we can define the thermodynamic force 
Xi = F/Tc and its conjugate thermodynamic ffux Ji = x 
for steady-state heat engines as well as Xi = —W/Tc and 
Jl = 1/Tcyc for cyclic heat engines. We can also define 
the other thermodynamic force X2 = 1/Tc — 1/Th and its 
conjugate thermodynamic flux J2 = Qh for both types of 
heat engines. By using these thermodynamic fluxes and 
forces, the power P is rewritten as P = —JiXiTc- We 
assume that these thermodynamic fluxes and forces satisfy 
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Fig. 1: Schematic illustration of the minimally nonlinear irre- 
versible heat engine described by eqs. (fTl)) and ((TJJ . 



the following extended Onsager relations: 

Ji = LiiXi + L12X2, 

■h = L21X1 + L22X2 — 7h. Ji^ 



(11) 
(12) 



where the features of the Onsager coefficients L^j's in 
eqs. ([3]) and (|4]) are assumed to hold also in eqs. (|TT|) and 
P^ . In eq. p^ . the new nonlinear term — is in- 
troduced into the standard Onsager relation eq. @ and 
7/1 is assumed to be a positive constant as 7/1 > 0. We 
also assume that no other higher nonlinear terms arise in 
eqs. pT|) and p2|) by considering that the coefficients of 
higher nonlinear terms are too small to have effective con- 
tributions to cr compared to —^hJi^, as will be seen in 
eq. P5|) . In our model, we also assume that Xi and X2 
are not restricted to small values, unlike in the linear irre- 
versible heat engine described by eqs. (jSj and (|4]), where 
the limit oi Xi and X2 —J- is taken. 

Here we introduce the heat flux from the cold heat reser- 
voir 



P-Qh = -JiXiT,-J2 = J3 



(13) 



Then we can rewrite the heat fluxes eqs. and into 
more suggestive and symmetrical forms as (see fig. [1} 

J2 = ^Ji + L22{l-q^)X2--/hJi^, (14) 



J3 



L2iTc 
LiiTh 



Jl-i22(l-9')^2-7cJl', (15) 



where q is defined in eq. ([7]) and 7c is defined as 



ill 



Ih- 



(16) 



Wc assume 7c to be a positive constant as 7c > 0. We 
can use eqs. P^ and P^ to describe the heat engines 
instead of eqs. pT|) and p2|). regarding Ji as the con- 
trol parameter of the heat engines instead of Xi since Ji 
and Xi are uniquely related through eq. ()lip when X2 
is fixed. We call the heat engines described by eqs. p4|) 
and p5|) (or the extended Onsager relations eqs. pT|) and 
P2p ) the minimally nonlinear irreversible heat engines. 
The term "minimally" implies that we take into account 



only— 7/1 Ji and — 7c Ji as the nonlinear terms. As we ex- 
plain below, those terms will turn out to be the inevitable 
power dissipations accompanied by the finite-time motion 
of the heat engines. 

The power P = —JiXiTc = J2 + J3 can be rewritten as 



L21 



Tc 



(17) 



by adding eqs. p4p and p5|) . Here we immediately no- 
tice that the second terms in eqs. P^ and psp do not 
contribute to eq. p7)) at all. They mean just the direct 
heat transfer from the hot heat reservoir to the cold one, 
which arises in the case of the non-tight coupling con- 
dition \q\ 7^ 1. Note that similar decomposition of the 
power like eq. p7|) is also given in 02]. The first term 
in eq. p7|) is proportional to AT through r/c", meaning 
the power generation due to the temperature difference. 
On the other hand, the second term in eq. p7p remains 
nonzero even when AT = and means the power dissipa- 
tion which should necessarily be consumed once the heat 
engine moves at a finite rate ( Ji 7^ 0) regardless of how 
small Ji is. The power dissipation results in the increase 
of the total internal energy of the heat reservoirs. When 
AT = 0, this is nothing but the effect of Joule heating, 
which is seen if we can consider that Tc/Tn and Ji in 
eq. p7|) correspond to resistance and electric current, re- 
spectively. 

To clarify the physical meaning of each term in eqs. p4p 
and p5[) in more detail, we rewrite the total entropy pro- 
duction rate cr = -Qh/Th - Qc/Tc = -.hl^h - -hlTc 
as 



a ^ L22{\ - q^)X2 



7 

- 7/> Ji'^2, (18) 

^11 



by using eqs. p4)) and p5|) . The first term means the en- 
tropy increase rate of the heat reservoirs due to the direct 
heat transfer. The second term comes from the inevitable 
work consumption due to the finite-time operation. The 
third term in eq. psp arises due to the presence of — 7/1 Ji^ 
in eq. p2p . In the case of the linear irreversible heat engine 
described by eqs. ([3]) and Q, the third term is suppressed 
and the non-negativity of g restricts the Onsager coeffi- 
cients iy 's to eq. ([8]). Even in our model, we assume that 
the restriction eq. (|S]) still holds although X\ and X2 are 
not restricted to small values. Then the non-negativity 
of eq. (jlSp is always guaranteed since it is rewritten as 
a = T22(l-<?')^2' + (7/./r/, + 7c/rc) Ji" by using eq. (HH): 
the first term is always non-negative due to eq. and the 
second one is also always non-negative due to the assump- 
tions 7/i > and 7c > 0. 

Efficiency at maximum power. We consider the 
efficiency of the heat engine 77 — W/Qh = P/Qh = P/J2, 
where P and J2 are given in eqs. p7)) and p^ . respec- 
tively. When X2 and 's are given, the maximum power 
is realized at Ji = L12X2/2 as a solution of dP/dJi = 
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and then we obtain the maximum power P* and the effi- 
ciency at the maximum power if as 



Vc 



2 2-q^l + r^c/mi+Jc/lh))) 



(19) 
, (20) 



respectively. The formula eq. (|20p is the main result of 
this paper. We notice that it includes the formula eq. ^ 
of the linear irreversible heat engine as the linear term of 
77c =i AT/T in the limit of AT 0. We also notice that 
eq. ([20|) has the lower bound 77' and the upper bound 
at a fixed q as 



q _ Vc q . * . '7c 



2 2-g2 



2 2 - g2 (1 + 7/c/2) 



Vl, (21) 



by taking the asymmetrical dissipation limits "fc/^h — >■ 
00 and jc/jh — >■ 0, respectively. Moreover 77^ takes the 
maximum value 



V- 



< 



Vc 



V+, 



(22) 



when the tight-coupling condition |g| = 1 is satisfied in 
rj'^. Therefore this 77+ is the upper bound of rj* for the 
minimally nonlinear irreversible heat engines. 

We note that 77+ has also been found in the previous 
studies on the efficiency at the maximum power: in a 
finite-time heat engine model [40], where the heat fluxes 
are assumed to obey a specific conduction law, 77_|_ arises 
as a limiting case. In a Feynman ratchet model [41] . 77_|_ 
has been obtained as Feynman efficiency under the no heat 
leak condition between the heat reservoirs. In a few finite- 
time Carnot cycle models [24l|34l[35] , 77+ has been found 
as the upper bound of 77* in the asymmetrical dissipation 
limit. Finally in [35] (see also [35]), 77+ has been proved to 
be the upper bound of 77* in the asymmetrical dissipation 
limit in a general and mo del- independent way. But the 
proof in [35] is limited to the case of stochastic steady- 
state heat engines. Our theory can be applicable to cyclic 
heat engines as well and unify these previous results. 

Here we stress physical importance of 77+: if the effi- 
ciency at the maximum power of a finite-time heat engine 
exceeds 77+, it implies that the heat engine works under 
higher degree of nonequilibrium. In fact, we can see that 
77* of the finite-time Carnot cycle model of ideal gas re- 
ported in [551155] exceeds 77+ due to the higher nonequilib- 
rium effect [15]. Therefore 77+ could be a criterion for de- 
termining the degree of nonequilibrium of finite-time heat 
engines. 

Example: low-dissipation Carnot engine. For 

a demonstration of the validity of our theory, we show that 
the low-dissipation Carnot engine [3S] is described by the 
extended Onsager relations eqs. PT]) and (|12p . Here the 
low-dissipation Carnot engine is a heat engine model pro- 
posed as a finite-time extension of the quasistatic Carnot 



cycle. It assumes the specific form of the heats transferred 
from the heat reservoirs during the isothermal processes 
as 



ThAS - 
-TcAS 



Tc^c 



(23) 
(24) 



where AS is the quasistatic entropy change inside the heat 
engine during the isothermal process in contact with the 
hot heat reservoir, t/j and Tc are the durations during the 
isothermal processes in contact with the hot heat reser- 
voir and the cold one, respectively, and Hh and Ec are 
positive constants. We consider that the constants 
and Ec contain the details how the engine deviates from 
the quasistatic limit. The assumption eqs. (|23p and p4| 
means that the lowest deviation from the quasistatic heat 
should be proportional to the inverse of the duration. In 
a stochastic finite-time Carnot cycle model [53] analyzed 
by the Fokker-Planck equation, such a term indeed 
arises [cf. eq. (16) in [24]]. In a finite-time Carnot cy- 
cle model of ideal gas analyzed by the molecular kinetic 
theory [5S], it is also confirmed that the lowest deviation 
is proportional to r^^ [cf. eq. (11) in [5S]]. Finally such 
a term also arises in a quantum dot Carnot engine 
model [34] analyzed by the master equation approach [cf. 
eq. (28) in [Ml]. Therefore the assumption of the specific 
form of the heats eqs. ([23]) and (f24| has microscopically 
been justified in these models. Additionally, we neglect 
higher order terms such as 0(t~^) in eqs. ([53]) and ((5^ 
in this low-dissipation approximation. 

We can express the power P = [Qh + Qc)/{Th + Tc) of 
this engine as P = (ATA5-r^S,,/T,,-T,I],/T,)/(r,, + r,) 
by using eqs. ([23| and ()24p . Maximizing this power by the 
durations Th and Tc as dP/drji = OP/Otc = 0, we find the 
physically relevant solutions as 



Th = 



^Th^h 



{Th - T,)AS 
{Th - T,)AS 




Th\ (25) 



(26) 



Then by using the definition 77 — {Qh + Qc)/Qh and 
eqs. ([231), ([24]), ([25]) and ([26]), we can obtain the efficiency 
at the maximum power 77* as 



(1 + ) 



(27) 



We can easily notice that eq. ((57)) is bounded from the 
lower side and the upper side as 



Vc ^ ^ ^ Vc 
2 - ' - 2 - 77c 



(28) 



by taking the asymmetrical dissipation limits 'Sc/'Sh — > 00 
and Sc/E^ — >■ 0, respectively [3S] (see also [33] for the 
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derivation of these bounds in a quantum dot Carnot en- 
gine model). In [35], it is stated that observed efficien- 
cies of various actual power plants tend to locate between 
these two bounds. It is also interesting to see that the 
same bounds were derived in a different finite-time heat 
engine model based on a specific heat conduction law [3D] . 
By comparing eq. (j28p with eq. (j2I[) . we may consider 
that the tight-coupling condition jg] = 1 holds in this low- 
dissipation Carnot engine. We can prove it by writing 
the extended Onsagcr relations of this engine explicitly as 
follows. 

First, we consider the total entropy production rate 



(T = — 



Qh Qc _ 

W 

'T4a + l)Th 



W 



1 



(29) 



where we have defined the parameter a as a = Tc/rh and 
the dot denotes the quantity divided by the one-cycle pe- 
riod Tcyc = Tfi + Tc = {a + l)Th. From the decomposition 
a = J1X1 + J2X2, we can define the thermodynamic forces 
Xi = -W/Tc, X2 = 1/Tc - 1/Th and their corresponding 
thermodynamic fluxes Ji = l/((a4-l)T/i), J2 = Qh- Using 
eqs. ([23]) . ([24]) and the definitions of the thermodynamic 
forces and fluxes, we can easily calculate the Onsagcr co- 
efficients Lij 's and the constant 7^, of this low-dissipation 
Carnot engine as 



ill = 

L12 = 
L21 = 
L22 = 



Tc 




+ TcEc/a)ia- 


hi)' 




ThTcAS 






+ TcT.c/a)(a - 


hi)' 




ThTcAS 






+ TcT.c/a){a - 


hi)' 




n^TcAS^ 





(30) 
(31) 
(32) 
(33) 



(r,,S,,+T,E,/a)(a + l)' 
jh^nJ:h(a + l), (34) 



respectively. 7c is also given by 

Tc^c{a + 1) 



7c 



a 



(35) 



by using eqs. ([201), ([H]) and ([H]). We notice that the 
reciprocity L12 — £21 surely holds from eqs. ([3T]) and 
(j32[) . Moreover we can confirm the tight-coupling con- 
dition \q\ ^ \Li2/VLiiL22\ = 1 from eqs. ^ and 
([33]) as expected. We can obtain 77* of the low-dissipation 
Carnot engine as 



Vc 



2- 1 + 2 1 



aTh'^h 



, (36) 



by substituting \q\ = 1, eqs. ([34| and ([35]) into eq. ([20 
We can also obtain P* as 



AS^AT^ 



MTh^h+Tc^c/a){a + l)-' 



(37) 



by substituting \q\ — 1, eq. ([55]) into eq. ([TO]) . However we 
notice that eq. ([57]) still contains the tunable parameter a 
and can further be maximized as dP* [a) / da = 0, which 
reduces to a = ^TcYjc/ {TjJ^h) = a* . From eqs. ([25]) and 
([26]) . we can see that a* = Tc* /rt* holds. Substituting this 
a* into eq. ([36| . we finally reproduce eq. ([27| . Therefore 
we can conclude that the low-dissipation Carnot engine is 
exactly described by the extended Onsager relations. In 
other words, the inclusion of the power dissipation term 
—IhJi^ into the Onsager relation as in eq. ()12p is justi- 
fied by this explicit example, whose assumptions eqs. (|23p 
and ([^^ are consistent with the microscopically analyzed 
models iMKMj. 



Summary and discussion. — We proposed the min- 
imally nonlinear irreversible heat engine described by the 
extended Onsager relations, where a new nonlinear term 
meaning the power dissipation is added to the heat flux 
from the hot heat reservoir in the standard Onsager rela- 
tion and no other nonlinear terms are assumed to arise. 
Thus our model can be regarded as a natural and mini- 
mal extension of the linear irreversible heat engine. We 
formulated the efficiency at the maximum power 77* of 
our model and showed that it is bounded from the up- 
per side by 77c/(2 — 7?c)- This upper bound can be at- 
tained when the heat engine satisfies the tight-coupling 
condition |g| = 1 and the asymmetrical dissipation limit 
IcMh — > is taken. As a demonstration of the validity 
of our theory, we explicitly wrote down the extended On- 
sager relations of the low-dissipation Carnot engine j35j 
and confirmed that it satisfies the tight-coupling condi- 
tion \q\ =■ 1. Though the low-dissipation Carnot engine is 
an example of the cyclic heat engine, we should note that 
the power dissipation terms arise also in a few steady-state 
systems [43]]45] , where analytical calculations of the On- 
sager coefficients 's, 7/1 and 7c are explicitly done based 
on a molecular kinetic theory. These calculations and the 
present example of the low-dissipation Carnot engine in 
this paper could support the validity of our theory, which 
treats the cyclic heat engines and the steady-state ones in 
the unified manner. It will be a future challenge to find 
the upper bound of the efficiency at the maximum power 
for more general irreversible heat engines with higher non- 
linear terms beyond our model. 
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